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Abstract 

This work concerns the loss of energy of a material system due to gravitational radiation in 
Einstein-aether theory — an alternative theory of gravity in which the metric couples to a dynamical, 
timelike, unit-norm vector field. Derived to lowest post-Newtonian order are waveforms for the 
metric and vector fields far from a nearly Newtonian system and the rate of energy radiated by the 
system. The expressions depend on the quadrupole moment of the source, as in standard general 
relativity, but also contain monopolar and dipolar terms. There exists a one-parameter family of 
Einstein-aether theories for which only the quadrupolar contribution is present, and for which the 
expression for the damping rate is identical to that of general relativity to the order worked to 
here. This family cannot yet be declared observationally viable, since effects due to the strong 
internal fields of bodies in the actual systems used to test the damping rate are not included. 
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I. INTRODUCTION 



That the world exhibits exact Lorentz invariance is a central hypothesis of modern physics, 
and as such, warrants scrutiny. Interest in testing this hypothesis has been growing in recent 
years. This is likely due to hints that this symmetry may be broken in various candidates 
for theories of quantum gravity and models of physics beyond the Planck scale, including 
strings jl|, loop quantum gravity Q], and noncommutative field theories [sj]. The review 4] 
discusses various theoretical models that feature Lorentz-symmetry violating effects and 
observational searches for violations. 

If we wish to incorporate Lorentz violation into a gravitational setting, we must employ a 
mechanism that breaks this symmetry while preserving the distinct symmetry of diffeomor- 
phism invariance. Einstein-aether theory is a variation of general relativity (GR) that does 
just this. Einstein-aether theory, or 'ae-theory' for short, is a classical, vector-tensor theory 
of gravity in which the vector field 'aether' is constrained to be everywhere timelike and of 
fixed norm. The aether can be thought of as an effective remnant of unknown, Planck scale, 
Lorentz-violating physics. It defines a "preferred" frame, while its status as a dynamical 
field preserves diffeomorphism invariance. The condition on the vector norm (which can 
always be scaled to unity) ensures that the aether just picks out a preferred direction and 
removes instabilities in the unconstrained theory (see J5|). A review of properties of this 
theory and references to earlier work, can be found in [6J. 

Much work on ae-theory has focused on deriving observational constraints on the theory's 
free parameters, denoted here 1,...,4). Recent efforts have determined the 

form of the theory's post-Newtonian expansion and the nature of linear, source free wave 
phenomena, and have demonstrated that the theory remains observationally healthy for a 
large range of c n values. The parameterized post-Newtonian (PPN) expansion jl,^ describes 
weak field, slow motion gravitational effects in terms of ten parameters that depend on the 
theory's free parameters. The PPN parameters have been constrained by solar system based 
tests to have the values predicted by GR, up to small errors. In ae-theory, only two of the 
ten parameters, «i and a 2 , differ from the GR values 0, [lo|, 11]. These can be set to their 



GR values of zero by special choice of two of the four c n . The propagation of plane waves 



on a flat metric, constant aether background was examined in 



121 ] . It was shown that there 



generally exist five independent wave modes — two spin-2, two spin-1, and one spin-0 — that 
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travel at three different speeds that depend on the c n . The energy densities associated with 



these modes were derived in [13J] and also depend on the c n . 

Constraints on the c n are then imposed by the requirements that a% and «2 vanish, 
and that wave modes have positive energy densities and real speeds. Another constraint 
arises from imposing superluminal speeds, which is required to a first approximation by 
the absence of vacuum Cerenkov radiation of gravity-aether shock waves [5(. Yet another 



constraint can be derived from primordial nucleosynthesis 



14j . It was shown in that all 



of these constraints are satisfied by a large two parameter family of ae-theories. 



Additiona 
begun in 



tests of ae-theory can come from the study of stellar solutions and black holes, 



l^ l , and from the "ultimate" test 



?| provided by observations of binary pulsar 



systems. This work begins that examination with a calculation of the generation of gravity- 
aether radiation by a nearly Newtonian source and the subsequent energy loss, or radiation 
damping, of the source. A formula is derived for the rate of change of energy: 

where Qij is the trace-free quadrupole moment of the source, / is the trace of the second 
moment, £j is a dipolar quantity defined below, and A, B, and C are dimensionless combi- 
nations of the c n ; Gn is the value of Newton's constant that one would measure far from 
an external gravitating source; the angular brackets indicate a time average over a period of 
the system's motion. This formula generalizes the "quadrupole" formula of standard general 
relativity, which predicts a similar expression but with A — 1, B — C — 0. 

In the case of a system of two compact bodies, this expression takes the form 

d£ //GV/im\ 2 /8 / 2 /^ r ^2^ ak>1^ T 



dt 



^((T)0-<« + <) ,+ «'))' < 2 > 



where /i is the reduced mass of the system, m the total mass, v the relative velocity of the 
bodies, and r their orbital separation, which I assume is much larger than the size d of the 
bodies; T> is the difference in self-gravitational binding energy per unit mass of the bodies, 
and the coefficient C is another dimensionless combination of the c n . 

This expression gives the lowest-order effects in a post-Newtonian (PN) (weak field, slow 
motion) expansion. Aside from the (GAr/xm/r 2 ) 2 prefactor, the first two terms are 0(GNm/r) 
and the last is 0{{G^m/ d) 2 ). It does not take into account strong field effects that may be 
important when the fields are not weak inside a given body. The strength of the field of a 



compact body can be characterized by the quantity {G^m/d)] this is "small" for the sun 
(~ 10~ 6 ) or a typical white dwarf (~ 10~ 3 ), and "large" for a typical neutron star (~ 10 _1 ) 
or black hole (~ 1). Thus, the field should be strong within the systems actually used to 
measure the damping rate. Strong field effects on the damping rate of a compact body can 
be associated with a dependence of the body's gravitating mass on the ambient non-metric 
fields — that is, with a violation of the strong equivalence principle jg]. These effects are not 
present in GR at lower post-Newtonian orders. Their presence in ae-theory will be examined 
in future work 17]. 



The damping rate can be tested by observing the rate of change of the orbital period 
P of various binary systems, since (dP/dt)/P = —(3/2)(d£/dt)/£, equating the energy 
radiated to minus the change in mechanical energy of the system. In practice, this test 
is conjoined with tests of other "post-Keplerian" (PK) parameters s|, in particular 
the rate of advance of periastron (the point at which the two objects are closest to each 
other) and the redshift or time delay due to the gravitational field of the system. These 
"quasi-static" [3] parameters are determined by the post-Newtonian forms of the fields 
and the effective equations of motion for the compact bodies. The conjoint technique is 
necessary, because the expressions for the PK parameters depend on the unknown masses 
of the systems' bodies. The expressions for the parameters will depend on the two masses, 
other measurable parameters, and a given theory's free parameters. Measurement of three 
mass dependent parameters, for fixed values of the theory parameters, gives three bands with 
widths due to errors in the two-dimensional space of mass values. The theory is consistent 
for those values of the free parameters if the bands overlap. The predictions of GR have 
been validated in this way using data from various binary systems containing pulsars, whose 



181 ] for details. 



regular pulsing provides an accurate measuring device; see the review 

For ae-theory, I find that if one assumes the strong field effects are negligible so that the 
results below are adequate, then there exists a one-parameter family of theories that satisfy 



all of the constraints summarized in 



111 ], and whose predictions for the PK parameters 



match those of GR to the order worked to here. This can be seen as follows. To lowest 
PN order and neglecting strong field effects, the quasi-static parameters can be determined 
within the PPN framework [7j. Consequently, when «i and a<i are set to zero, so that 
all of the ae-theory PPN parameters match those of GR, the two theories will make the 
same predictions for the quasi-static parameters. In this case, the c n can be constrained by 
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requiring that the damping rate equal that of GR. But as shown below, the radiated fields 
contain only quadrupolar contributions when a\ and «2 are set to zero. The damping rates 
then coincide when A is set to one, which can be done by imposing one condition on the 



two remaining free c n . To be consistent with the observational tests summarized in [11], this 
curve of theories must intersect with the allowed two-parameter family demarcated there 
(and below in Sec. |V]). This is the case all along the curve, as long as c_, c + > 0. 

The calculation will proceed as follows. In Sec. [Til a weak field expansion of the field 
equations is performed. The perturbations are shown to satisfy the wave equation, with 
matter terms and nonlinear terms acting as sources. In Sec. IIH[ these equations are solved 
via integration of the sources with Green's functions. The source integrals are approximated 
in terms of time derivatives of moments of the sources, and evaluated to order of interest 
using the PPN expansion of the fields. In Sec. IIV| an expression for the rate of change of 
energy contained within a volume of space is defined, and evaluated in terms of the wave 
forms. I will conclude with a discussion of the constraints on the c n implied by observations 
from binary pulsars, in Sec. IVl 

I use the conventions of [l9[. In particular, the metric signature is (—,+,+,-1-). This 
differs from the published version of this article 2(|; however, this signature is vastly more 
convenient for performing calculations involving a time-space decomposition, and should 
have been used from the start. The coefficients in the action below will be defined so that 
the results for the PPN parameters and the damping rate match the published results. 
Spatial indices will be indicated by lowercase Latin letters from the middle of the alphabet: 
k, . . . . One exception is that the coefficients 01,2,3,4 will often be referred to collectively 
as c n in the text, when no confusion should arise. Upon performing the weak-field expansion 
of the field equations (Eqn. ffl3|) onwards), indices will be raised and lowered with the flat 
metric r] ab . Repeated spatial indices will be summed over, regardless of vertical position: 
Tu = X/i=i 3^&- The flat-space Laplacian will be denoted by A: Af = f\a. Time indices 
will be indicated by a 0; time derivatives will be denoted by an overdot: / = dof. I adopt 
units in which the flat-space speed-of-light c = 1. 
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II. FIELD EQUATIONS 

In this section, I will expand the ae-theory field equations about a flat metric, constant 
aether background, obtaining a set of wave equations with matter terms and nonlinear terms 
as sources. 



A. Exact equations 

I begin with the standard four-parameter ae-theory action S 

S = J d 4 x y/\g\ (R - K ab cd V a u c V b u d + Ku a u b g ab + 1)) , (3) 

where 

K a \ d = {c l9 ab g cd + c 2 6 a c 6 b + c 3 6 a d 5 b - c,u a u b g cd ) . (4) 

The sign of the c 4 term looks awkward, but it allows for easy comparison between results in 
this version of the article and in the published version {^(J. I will use shorthand for certain 
combinations of the c n : 

C M = Ci - C 4 , C123 = Ci + C 2 + C 3 , C± = Ci ± c 3 . (5) 

In addition to the ae-theory action, there is an aether-independent matter action. The 
matter can be assumed to couple universally to some metric since Lorentz-violating effects 
in nongravitational interactions are already highly constrained {4, 21]. Aether couplings are 
then excluded from the matter action, and the field g a b is identified as this universal metric. 
The resulting equations of motion consist of the Einstein equations 



where 



with 



Gab — Sab — S^GTab, (6) 

G a b = Rab — -Rg a b, (7) 

Sab =V c (K c {a u b ) + K(a h) u c - K (a c u b )) 

+ ci ( V a u c V b u c - V c u a V c u b ) + c 4 {u c V c u a ) {u d V d u b ) (8) 
+ Xu a u b - ^g ab (K c d V c u d ), 

K\ = K ab cd V b u\ (9) 



and T a b is the matter stress tensor. There are also the aether field equations, 

V b K\ = -\u a - c A {u c V c u b )V a u\ (10) 

and the constraint 

g ab u a u b = -1. (11) 
Eqn. (flQjl can be used to eliminate A, giving 

A = u c V a K\ + Ci (u c V c u a )(u d V d u a ). (12) 
B. Linear-order variables 

I will now expand the exact equations about a flat background. I assume a Minkowskian 
coordinate system and basis with respect to which, at zeroth order, the metric is the 
Minkowski metric r] ab and the aether is purely timelike. I then define variables h ab and 
w a , with 

h ab = gab ~ Vab, W° = U° - 1, W % = u\ (13) 

I assume that h ab and w a fall off at spatial infinity like 1/r. 

I will further define variables by decomposing the above into irreducible transverse, or 
"divergence- free" , and longitudinal, or "curl-free", pieces. The decomposition is unique 
and well-defined in Euclidean space, having imposed the above boundary conditions (one 
is essentially solving Laplace's equation — see 22j for more discussion). First, consider the 
spatial vectors w l and h^i, and define the following variables: 

hot = 7i + l,i = v i + v> i , (14) 

with j = v\ = 0. Next, consider the spatial components of the metric hij. A symmetric, 
2-index tensor on Euclidean space can be uniquely decomposed into a transverse-trace- 
free tensor, a transverse vector, and two scalar quantities representing the transverse and 
longitudinal traces: 

fkj = (pij + 7^Pij[f] + 2<P(i,j) + <f>,a, (15) 

where 

= 4>%j,j = 4>jj = 4>i,h (16) 
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and 

I';;!. 5 {j Af - A,: (17) 
hence, Pij[f]j = 0, and ha = A(f + <f>). Further define 

F = Af. (18) 

The list of variables then consists of a transverse-traceless spin-2 tensor faj, transverse spin-1 
vectors ji, u l , fa, and spin-0 scalars h 00 , w°, 7, v, F, (p. 

I will impose coordinate gauge conditions below, after expressing the field equations in 
unfixed form. The standard gauge to impose when performing the analogous calculation in 
conventional GR is the "harmonic" or "Lorentz" gauge, 2h ab ,b = if d h c ^ a , as this happens to 
reduce the field equations to a simple form when they are expressed covariantly. Some variant 
of this condition has a similar effect in several other alternative theories of gravity, as seen 



in [231 ] . Here, the increased complexity of the equations and the noncovariant decomposition 
of them and the field variables means that no obvious extension of the harmonic gauge has 
such a utility. Instead, the gauge will be chosen somewhat arbitrarily so as to eliminate 
certain variables: 

= = h 0i ^ = hifak]h (19) 

or equivalently, 

= ^ = 7 = 0,. (20) 
An infinitesimal coordinate gauge transformation has the linear-order form 

tf/i«6 = £«,& + 6,« 5w a = -?. (21) 

The conditions ( 1201) can be realized while in an arbitrary gauge (a prime denotes that the 
variables are evaluated in the original gauge) by choosing £ = — (7' + v') and the transverse 
part of £j as — and by solving for the longitudinal part £ of via £ = v' . One constraint 
on the choice of gauge is that it must be a valid PPN gauge, as defined in so that the 
integrals of Sec. IIII CI can be evaluated by expressing the variables in terms of their PPN 
expansion. The above is a valid, albeit nonstandard, PPN gauge (in contrast to the gauge 



chosen in 



12 



I*- 
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C. Linearized equations 



I now express the field equations ([6]) and (jTOj) in terms of the above variables, and arrange 
them in the form 

Gab - S ab = 8nG(T a b + t ah ) , (22) 
K\ b = 8nGa a , (23) 

where the overbar denotes the portion of the tensor linear in h a b and w a , and the nonlinear 
source terms t a b and a a are defined to a given order in the variables by asserting that the 
above equations equal the exact equations to that order. It will prove convenient to combine 
the equations in the form 

Gab - Sab + S[ a K c b ] tC = 8nGT ab , (24) 

thus defining the source 

Tab = T ab + tab + tf a (T b] . (25) 

Identities satisfied by the linear-order terms will imply conservation of r a b- 
The constraint (TTTj) to linear order is 

w° = h 00 . (26) 

I will use this result to eliminate to . The form of nonlinear terms will not be needed as 
explained in Sec. IIII CI 
Now, 

Gab = -^{^h ab ~ hab) ~ ^Kab + Kiatf + 7,Vab(&h -h~ h c / d ), (27) 

where h = r] ab h a b- Hence, 



Gij = -\ [A0y - 4-] + [0 W) - 7 (iJ) ] + \Pa\M -f- 2h 00 - 24> + 4 7 ] - - i J., r (28) 

S« = ~A( 74 - ^) - (29) 

Goo = ~AF. (30) 
The linear-order forms of the covariant derivatives of u a are 



Vow* = w l + hoi - -hao,i 

= v l + 7^ + (z> + 7 - ^^oo),* 



(31) 



1; 
- 
2 

-fai + v]i + l[ hi ] + 0o-,i) + -rPiAf] + iy + ^)«, 



and V a u° = 0. 
From 

follows 



1 



S w + -KS = KM d = c + d s (VW)) 



l&(c + (S + + ((c+ + c 2 )(2v + 0) + c 2 /) .). 



S 0i - -K ai ' a = K 0l + K MJ = c u d (V u l ) + c_a,-(VM) 

c_ 1 • 

= c i4 (i> 1 + 7i) - tt A (^ + Ti) + ci 4 (i> + 7 - o^oo) 
z z 



and 



(32) 



Sab = K(ab) + S (a K b) C ,d ( 33 ) 



5^ = d (c+V(V) + c 2 c%V fc w fc ) 

= ^ + c + (z>^ + + ^[c 2 (2z> + # + /) + ^/] (34) 
+ ^((c 2 + c + )(2z>+0)+c 2 /).., 



(35) 



(36) 



S 00 = Cudj(V ui) = a(c 14 (u + 7 - ^oo))- (37) 
The above expressions indicate that the linear-order terms satisfy the identity 

{G ab -S ab + 5^ a K c blc )' b = 0. (38) 

This implies that the source r ab ( 1251) obeys a conservation law 

r ab ' b = T ai>i - foo = 0. (39) 

D. Wave equations 

The above equations can now be decomposed as the variables. The field equations for 
variables of different spin will separate. I will impose below the gauge conditions f[2"0"j) . The 
following results are equivalent to those of [3] expressed in a different gauge when T ab = 0; 
in particular, the count of independent plane wave modes — two spin-2, two spin-1, and one 
spin-0 — and the expressions for the wave speeds are recovered. 
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I now consider the different spins in turn. Define the set of operators 
for % = 0,1,2. 



1. Spin- 2 



The transverse-traceless part of the space-space components of ( 124T) gives 

□ 2 o - = -167rGrJ T , 

with 

(* 2 ) 2 ' 



l-c + ' 

and where TT signifies the transverse-traceless projection. 
2. Spin-1 

Now the spin-1 variables. The transverse parts of f l24l) give 

A(c+i/ i + 7i ) = -167rG^, 

and 

c u (i? + 7i) - ^ A ( c -^ + C 1 - C -H) = -&kGt£- 
where the T signifies the transverse projection. These relations imply 

+ 7.) = ~ 167r ° (c + r i0 + (1 - c + K) T , 

LC\ — C_|_C_ 

with 

2 2ci — c + c_ 



2(1 - c + )c M 
5. Spin-0 

Now consider the spin-0 variables. The transverse-trace and longitudinal-trace port 
of the space-space components of fl24|) give 



[l + 2c 2 + c+)F -A(F- 2h 00 - 2(1 + c 2 )0) = -167rGr, 
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a i 



and 

(1 + c 2 )F + ci 23 A0 = -lQnGr^ (48) 
where rjj = Ty — r?j. The time-time component of (|24p gives 

A(F-c 14 /i o) = -167rGr 00 , (49) 

and the longitudinal space-time component of (124"]) gives 

A((l + c 2 )f + c 123 0) . = -16^, (50) 

where = r i0 — r^. These equations imply 

1QttGc 14 2 + 3c 2 + c+ l 2 . , , 

□ F = -— (r« 7- + — Too), 51 

2 - C M C123 Cu 

with 

(2 + 3c 2 + c + )(l - c + )ci4 
Further implied by these and the untraced, transverse-trace part of (I24"|) is the equation 

□o/,ii = tJ. (53) 

The form of the source r^- is unimportant; only the fact that satisfies a sourced wave 
equation is needed so that later eqn. (I9"8l can be applied when evaluating the damping rate 
expression in Sec. HV] 

III. EVALUATION OF SOURCE INTEGRALS 

The above equations can be formally solved via integration of the sources with the appro- 
priate Green's function, and the resulting integrals approximated in terms of time derivatives 
of moments of the source. Upon doing so, the nonstatic contributions to the fields to desired 
accuracy at points far from the material source depend on two integral quantities, the second 
mass moment of the material source 



/, = J top**,, (54) 
where p = T 00 to lowest order, and the integral 

Ei= / d 3 xai, (55) 



where <7j are the quadratic terms from the aether field equation (1231) . 
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A. Approximation of source integrals 



Equations of the form 

- (s)~ 2 ^ = -167TT, (56) 
can be solved with outward-going disturbances at infinity by writing 



3 r(t-z/s,x') 



if>(t,x) = 4j d A x'^ y \ (57) 

where z — |x — x'|. 

The source integral can be simplified with a standard approximation [7j] . As indicated by 
the energy loss rate expression ( 11011) . only the portion of the fields that fall off as (1/r) are 
of interest. A weak field, slow motion assumption will be made: the material source should 
be described by a mass m, a size L, and a time-scale T such that (Gniti/L) and (L/Ts) 
are small quantities. Then only terms of interesting order are retained in the following 
expansion: 

^ x )-^(E^^/^-^ x ')(^r)> (58) 

m=0 ' J 

where R = |x| and x = x/i?, and R~^> L. 

Now, the following sleight of hand justifies solving the decomposed ae-theory equations 
by first approximating the integral on the right side of (1571) using the full r and then taking 
the projection. Introduce the notation [[r(t, x')]] = r(t — z/s,x.'). Because the quantity on 
the right side of (15 7p depends on x only through z, it follows that 

d [Mt,x')]] _ f 8 ,[[r(t,x')]k , mV(t,x')]] 



dx l J z J dx' 

It follows from this that, e.g., 



, n[9/r(t,xQ]] _ (5g) 

Z J z 



[bg^oij _ n M^i \\ (60) 

after discarding integrals of total derivatives, where T on the left side signifies transverse with 
respect to x', and on the right side transverse with respect to x. As a further convenience, 
it follows that to 0(1/ R), the transverse projection is equal to the algebraic projection in 
the direction orthogonal to x. 

Additionally, there are the Poissonnian equations (H3|) . (j49p . and ([50]) of the form 

= -I671T. (61) 
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Solving via Green's function and expressing to 0(1/R) far from the source gives 




(62) 



The integrals of the sources in these particular equations happen to be conserved quantities. 
Thus, ignoring static terms in the wave forms, the equations are effectively imsourced, and 
thus imply 



to sufficient accuracy 

B. Sorting of source integrals 

To evaluate the integrals indicated by Eqn. fl58l) . I will express the sources in terms of 
their post-Newtonian expansions. I will further assume that the system is composed of 
compact bodies of individual size d L that exert negligible tidal forces on each other. 
The system will then have an orbital velocity v ~ y 'G^mjL. Following the discussion 
in the leading-order terms in the fields will be 0(G^m 2 /L), which give the quadrupolar 
and monopolar contributions, and 0((Gn17i 2 /d)v ), giving the dipolar contribution. Terms 
of these orders can only result from integrals of terms that are, respectively, 2PN and 2.5PN 
order. Integrals of interest can be identified by noting that since the rate of change of the 
system is governed by its velocity, assumed to be .5PN order, taking the time derivative of 
a quantity effectively multiplies it by a factor of v and raises it by .5PN orders. Also, only 
nonstatic, or non-conserved, terms are of interest as only the time derivatives of the fields 
will appear in the expression for the energy loss. 

I begin by considering the moments of r%j. First, the conservation law implies 



where the last two equalities hold to desired order and for the last I have used the Eulerian 
continuity equation for the fluid 



ip = 



(63) 




(64) 



(65) 



assumed to hold at 0(1.5). Then, 




(66) 
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which is of uninteresting order, as are remaining moments. 

I now consider the moments of Tjo- The integral of Tjo is conserved, so I ignore it. Next, 

TiDXj = - J n, = —iy. (67) 

The other moments are of uninteresting orders. 

I now consider moments of Too- First, the integral of Too is conserved, so I ignore it. Then, 



Too^i = - J m = (68) 
where the second equality ignores the static integral of To- Finally, 

J TooXiXj = Iij, (69) 

to desired order. 

C. Evaluating Sj 

I now consider the moments of o\. The terms in <j{ are at least 2.5PN order. The only 
integral of interest is thus Sj = J <Ji, which is 0({G^m 2 / d)v). At this point, I can explain 
why the nonlinear terms in the unit constraint ( TTTT) can be ignored. The previous subsection 
makes clear that only their appearance in cij need be considered. As follows from the PN 
forms given in Chapter 4, the nonlinear constraint terms are integer PN orders starting with 
2PN and have no free indices, and there are no field variables that are .5PN order. It follows 
that any nonlinear constraint terms appearing at 2.5PN order in G{ must do so in the form 
(terms) fli. Total derivatives do not contribute to so these terms can be ignored. 

I will evaluate Sj explicitly by expressing the fields in terms of the PPN expansion, but 
in the nonstandard coordinate gauge (1201 . The PPN forms in the standard PPN coordinate 
gauge, adjusted to the conventions used here, are reported in [ll] as 



<f>ij = (70) 

7, = - — S = -$±^.(y i + W i ), 0, = O (71) 
c_ 4 

h 00 = -A X , / = 20 = -2 X , (72) 

1 (R , v . 2ci + 3c 2 + c 3 + c 4 . 
7 = -7(6 + 01 -2a 2 )x, v= x, (73) 

4 ^Ci23 
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where 

Vi(x) = G N J d 3 x ,P ^, Wi(x) =G N J d 3 x' p{x)v f> Z \ (74) 
x (x) = -G N J d 3 x' p(x')z, (75) 
with z l = x l — x H (note that x — (K — W 7 *),*), and 

G " = T^tftY < 76 > 

2ci — c + c_ 

«i (ci + 2c 3 - c 4 )(2ci + 3c 2 + c 3 + c 4 ) , . 

" 2 = "o ro — "TX: • ( 78 ) 

2 (2-Ci 4 JCi23 

Adjustment from the standard to the nonstandard gauge is done by defining a gauge pa- 
rameter £ a with £ = — (7' + ^')> ^ — v \ii where 7', z/ are the standard-gauge values. Then 
in the nonstandard gauge, the variables are as above except v = 7 = and 

(c L +2^ ±£! ) 

Cl23 

With these forms, Ej can be evaluated, and after some algebra gives 

Xi = \jp((<x 1 -a 2 )V i + a 2 W$. (80) 

I will later consider the special cases of a single, compact, spherically symmetric body and of 
a pair of compact bodies that are static and spherically symmetric in their own rest frames. 
In the first case, spherical symmetry implies that Ej vanishes — otherwise it would define a 
symmetry-breaking spatial vector. In the second case, or more generally with n such bodies, 
Ej can be simplified via the following observations. First, the Newtonian potential U felt 
at a given body contains an 0{G^m/L) contribution from the presence of the other bodies, 
plus an 0{GNfn/ d) self-contribution U, 

U a (x a ) =G N [ d 3 x' P _ (81) 

J a l X a x I 

where the integral extends just over the "a-th" body. Second, spherical symmetry of each 
body implies that (£l a )ij — (l/3)^<Aj) where 

(n.)« = ~\g n [ Sxd 3 x' p{x)p{x ^ Z \ (82) 
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and 



1 r„ 



tt a = (tt a ) tt = -- I d s xpU a . (83) 



a 



Third, staticity of each body implies that Vi = ^2 a (v a yU, and similarly for Wi. These facts 
imply that 

/ pVi = 3 f P W t = -2 ^2(v a yn a , (84) 

J J a 

plus terms of 0(Gniti 2 v/ L). Therefore, to interesting order, 

^-(ai-I^X^U. ( 85 ) 



3 

a 



D. Wave forms 



I can now express the nonstatic, radiation-zone fields, to desired accuracy. For spin-2, 



where 



For spin-1, 



^' = ^fe) TT (86) 



Qij — hj — qI, I — la- (87) 



i -2G 

v l = 

R 2c 



— L ° + x QijXj - 2S,) T , (88) 

l - c+c_ V(l - c+)si J 

H = ~c+v\ (89) 



For spin-0, 



where 



F = ^7T^(3(Z - 1 )-r'Q, r ' J +ZI+ — S^), (90) 
it 2 — C14 V c 14 s ' 

ho = —F, (91) 

C14 

h = Jl± ^U 02) 

Cl23 



z= 2(°i-**)(l-^). (93) 
3(2c + -c 14 ) V ; 
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IV. ENERGY LOSS FORMULAS 



I now turn to the expression for the rate of change of energy contained within a volume 
of space. Such an expression can be derived via the Noether charge method for defining the 
total energy of an asymptotically-flat space-time [24j, using the ae-theory Noether current 
derived in j^j. One can equivalently work in terms of pseudotensors, using the results 



of 



13| 



Following the discussion in the appendix of [24} , an expression for the total energy £ 
contained in a volume of space, for a theory linearized about a flat background, is given by 
the integral over that volume of a certain differential 3-form J a b c = J d tdabc- J can be obtained 
from the quadratic order ae-theory Lagrangian modulo a boundary term. It will depend on 
the metric, aether, and an arbitrary background vector field. To define the energy, choose 
the background vector field as t a = (d/dt) a . Choose the volume V to be that contained 
within a sphere of coordinate radius R. Then, 



£ 



C t 3[t] 



J It] 



d 3 xJ°[t], 



d(t-J[t}) 



d£lR 2 x l J l [t] 



(94) 
(95) 



R 



where in the second line I have used the formula Ct J = d(t ■ J) + 1 ■ dJ and the fact that 
dJ = when the dynamical fields satisfy the equations of motion 241 ]. 

I will define J with respect to the ae-theory Lagrangian L modulo a total derivative: 

1 



a =£ 



16nG 
IQtvG 



The procedure of 24j gives: 

r 1 



[hc(h 



ab,c 



f W\{Kb9 ab - r b ab g ac ' 

9 ah {^ ad T d cb -T: d T d ab )-K\Vj 



\h bc ' a + \v bc (h d / - h ad d 



(96) 



Kir; (7 " v 2 ' 2 ' ' " ' a/ 2 

- u a K bc - 2K [ab] u c ) - 2u b K a b ] - t a C 



d.b ca 

9 d V 



(97) 



where all indices on h a b are raised with the flat metric r] ab . 

I will presume that only the time average of the damping rate need be determined. It is 
then crucial to note that the damping rate is calculated to lowest nonvanishing PN order by 
treating the system as exactly Newtonian. The motion of the system can then be decomposed 
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into a uniform translation of the center of mass-recall the conservation of J Tjo-and a fixed 
Keplerian orbit in the center-of-mass frame. As indicated below, the field wave forms do 
not depend on the center-of-mass motion. It then follows that if the time average is taken 
over an orbital period, total time derivatives in (j9"7|) do not contribute. 

It is also useful to note that approximation of the source integrals (1581) implies that to 
0(1/R), 

jfjj = -(l/s)ipx j , (98) 
for field ij} satisfying the sourced wave equation (1561) . This relation then implies that 

= x%(x) = x%{x) = x'Pijifix)}. (99) 

These facts permit manipulation of terms within the integral, e.g.: 

J < h jk (f)j k i >= J < -y<fijk<fi,jiXk >= J < ~^A0A0Xj >, (100) 

where the angular brackets denote the time average. 
The energy loss rate then evaluates to 

£ = -L / dOietl^ + (2C1 " C+C - )(1 - C+ W + ^FF). (101) 
167rG J r \2s 2 si 4ci 4 s / 

The sign of the coefficient of the term for each spin is opposite to the sign of the energy 
density associated with linearized plane waves, as found in [131 ] . Thus, a positive energy 
mode implies energy loss due to radiation of that mode. 

The energy loss rate can be further evaluated by substituting in the expressions for the 
fields given in Sec. HHP} and performing the angular integral, using the results 

i- f<m^ = ^ (102) 



h I = ^ (2c-U)' (K(T^k)'(W + IM)- < 103 > 

substituted into expression (110 II) gives 

5 



£ = —G N {-[QijQiA+Byi I)+C[ SjE, ) ) , (105) 
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where 

/ c 14 wl 2cmc| 1 3(Z-l) 2 c 14 1^ . 

^ = I 1 " TJ U + (2c 1 -c + c_)^ + 2(2 -c 14 ) ^ (1 ° 6) 

B = ^f, (107) 
C = ^(^Al (108) 

0C44 V "^0 

where Z is given in (1931) . This constitutes the generalization to ae-theory of the quadrupole 
formula of general relativity, and contains additional contributions from monopolar and 
dipolar sources. 

The presence of the monopolar term means that a spherically symmetric source, such as 
a spherically pulsating star, can radiate at this lowest nontrivial PN order in the presence of 
the aether, whereas it would not in pure GR. In this case, Iij = (l/3)<5yJ, and as observed 
above, = 0. Only the spin-0 radiation fields are nonvanishing, and the energy loss rate is 
£ = —GnB(^(I) 2 ^. Bounds discussed in Sec. |Vl however, require that the PPN parameters 
«i ( !77|) and ot 2 ( I78|l . hence B (I107p . vanish for observationally viable ae-theories. 

For a binary system, treating the two bodies as static and spherically symmetric in their 
own rest frames leads to 



■(4f(iVj) — 3fifjr), (109) 



where m = mi + m 2 , \i = mim 2 /m, r = ri — r 2 is the relative separation and v = r. Then, 

^A = ^(^)W-llr 2 ). (HO) 

Also, 



= (ai - ^a 2 )—^—'Dxi, (111) 



2 Giyfim 

3 r 

where T> is the difference in binding energy per unit rest mass: 



v= Sh_£h (n2) 
m a m b 



Therefore, 

/ G N jjim\ 2 ( 8 , 2 \%\ A nf-\ 2 1 ^ 



£ = ~G N ((^p) (^(12^ 2 - H(r) 2 ) + 4£(r) 2 + (a, - ^ 2 ) 2 CP 2 )}. (H3) 
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V. PARAMETER CONSTRAINTS 



I will now discuss bounds on the c n that can be derived by imposing the observational 
constraints summarized in [11] and by comparing the damping rate prediction (j!13p with 
measurements of binary pulsar systems. The four c n can be reduced to one free parameter 
by requiring that the PPN parameters a\ and «2 vanish, and that the damping rate coincide 
with that of GR to lowest order. The theory will then satisfy all solar system based tests; it 
is not correct, however, to say that it would pass the binary pulsar test. This is because the 
fields inside a neutron star pulsar or black hole companion are not weak, and strong field 
corrections to the quasi-static parameters may arise. Nevertheless, the weak field results are 
adequate for small enough c n , as discussed in [17|. Therefore, it is useful to check whether 
this curve of ae-theories intersects the region allowed by positive energy, real frequency, 
vacuum Cerenkov, and nucleosynthesis constraints. 

The PPN parameters oti f!77j) and a 2 f!78l) for ae-theory were determined in llj. It was 
shown that they can be set to zero, so that all of the ae-theory PPN parameters coincide 
with those of GR, with the choices 

C2 = - 2cf+ ; xC3 ~ c l C4 = a. (H4) 

The positive energy, real frequency, vacuum Cerenkov, and nucleosynthesis constraints can 
then be satisfied if c\ and C3 lie within the region 

< c+ < 1, < c_ < — — — -. (115) 
+ 3(1 -c+) 1 ' 

When q;i and 02 vanish, so does Z (1931) hence B (11071) . and Sj. The fields then contain only 
a quadrupole contribution, and the ae-theory damping rate (11131) will match that of GR 
when A — 1. Solving numerically shows that a solution curve exists in (c + ,c_) space that 
intersects the allowed region (11151) for all positive values of c_; see Figure [TJ Thus, there 
exists a one-parameter family of ae-theories which satisfy all of the constraints summarized 
in HI], and which predict a damping rate identical in the weak field limit to that of GR. 

Observational error allows this curve to be widened into a band. As explained in the 
Introduction, the standard method of measuring radiation damping is to observe the rate of 



change of orbital period P of a binary system 



18], which will be proportional to S. The 



smallest relative observational uncertainty in P, which equates with the relative uncertainty 
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FIG. 1: Class of allowed ae-theories, if strong field effects in binary pulsar systems can be ignored. 
The four-dimensional c n space has been restricted to the (c+,c_) plane by setting the PPN pa- 
rameters a% and «2 to zero via the conditions (|114p . The shaded region is the region allowed by 
collected non-binary constraints, demarcated in (|115p . The dashed curve is the curve along which 
binary pulsar tests will be satisfied, assuming ae-theory weak field expressions. Along both this 
curve and the boundary of the allowed region, c_ — > oo as c+ — > 1. The curve remains within 
the allowed region for all c+ between and 1. Strong field effects may lead to system-dependent 
corrections to the binary pulsar curve for large c n ; however, all such curves will coincide with the 
weak field curve for \c n \ < (0.01) given current observational uncertainties 201 ] . 



in S, is of order 0.1% for the Hulse- Taylor binary B1913+16 [8lll8l|. This uncertainty permits 
the band \A — 1| < 10~ 3 . Numerical results indicate that at least for small c±, this band 
corresponds roughly to c± within about 10 -3 of the A = 1 curve. 
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